Sequential measurement-based quantum computing with memories 
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We introduce a general scheme for sequential one-way quantum computation where static systems 
with long-living quantum coherence (memories) interact with moving systems that may possess very 
short coherence times. Both the generation of the cluster state needed for the computation and its 
consumption by measurements are carried out simultaneously. As a consequence, effective clusters 
of one spatial dimension fewer than in the standard approach are sufficient for computation. In 
particular, universal computation requires only a one-dimensional array of memories. The scheme 
applies to discrete-variable systems of any dimension as well as to continuous-variable ones, and 
both are treated equivalently under the light of local complementation of graphs. In this way our 
formalism introduces a general framework that encompasses and generalizes in a unified manner some 
previous system-dependent proposals. The procedure is intrinsically well-suited for implementations 
with atom-photon interfaces. 



PACS numbers: 3.67.Lx, 03.67.Ac 

I. INTRODUCTION 

Arguably the most ambitious enterprise of quantum 
information science is the realization of a universal quan- 
tum computer. That is, a machine able to efficiently pro- 
cess quantum information in any desired fashion. Histor- 
ically, this quest was first pursued inspired by the clas- 
sical notion of an array of logical gates -a circuit- that 
sequentially act on a physical system encoding the logical 
state of the computer [T]. In this approach, the circuit 
model, all gates are decomposed into one-body (local) 
and two-body (entangling) unitary operations, and the 
output is read out by measuring the system's final state. 
A major advantage of the circuit model is that it requires 
the generation and coherent preservation of many-body 
entanglement only on a number of system constituents 
linear with the logical inputs. In what follows, we refer 
to this as Property PI. On the other hand, the sequence 
of operations must be adapted to each given computa- 
tion, which represents a considerable obstacle from an 
experimental viewpoint. In particular, entangling opera- 
tions must be actively applied on-line -during the course 
of the computation- while simultaneously preserving the 
system's quantum coherences. 

Recently, a conceptually different paradigm has been 
put forward: measurement-based quantum computation 
(MBQC) P]. There, the computation is performed by 
adaptive local measurements on certain quantum lat- 
tices in universal many-body-entangled resource states. 
The clear advantage is that no demanding on-line op- 
erations are required throughout the computation: only 
measurement bases must be (locally) adapted, whereas 
the generation of the resource state is independent of the 
given computation. More specifically, information pro- 
cessing proceeds via local measurements only (Property 
P2) ; and the architecture of the interactions needed for 



the resource-state generation is independent of the given 
computation (Property P3). The price to pay however is 
that universality is attained only with two-dimensional 
lattices, implying that the number of system compo- 
nents typically scales quadratically with the logical in- 
puts. To make things worse, such many-body entangled 
states seem not to be trivial to find in natural systems 
[3J. All in all, and despite the success in singling them 
out as ground states of specific, relatively simple Hamil- 
tonians [H [3] , their generation or coherent preservation 
still remain a challenge. 

In this work we introduce a hybrid model of quan- 
tum computation that combines the main practical ad- 
vantages of the circuit and MBQC models, and satisfies 
concurrently all the three properties P1—P3 listed above. 
The key idea of the construction is to apply the concept of 
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FIG. 1: Pictorial sketch of our sequential model for quan- 
tum computing. Long-lived quantum systems nij -quantum 
memories- iteratively interact at fixed interfaces with mov- 
ing registers -flying registers (curves). Registers do not need 
to be temporally synchronized and can also mediate the in- 
teraction between memories (thin curves). At each step, a 
graph state is generated and part of it immediately consumed 
by measurements at detectors Dj . Several physical platforms 
are envisaged for implementations in both discrete and con- 
tinuous variable systems (see text). 
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local complementation of graphs [SHE] to the most stud- 
ied MBQC architecture, namely the one-way quantum 
computer (1WQC) J2]. The resulting model preserves the 
advantages of 1WQC, while it relaxes the main obstacle 
toward its implementation. The procedure is completely 
general and can be applied both to discrete- variable (DV) 
systems of any dimension [31 [5] and continuous variables 

(CVs) pro]. 

Our proposal consists of an array of static, long-lived 
quantum registers -quantum memories that sequen- 
tially interact with moving, short-lived quantum registers 
-flying registers- via a fixed interaction (the only entan- 
gling operation required), see Fig. [I] At each step of the 
sequence, a graph state [7] (defined below) among mem- 
ories and flying registers is generated and, immediately 
afterwards, the flying sub-graph is consumed by adap- 
tive local measurements. The architecture of the system 
operations stays the same throughout and only the mea- 
surement bases must be adapted at each step, preserving 
the main advantages of the 1WQC. However, in contrast 
to the latter, here quantum information does not flow 
from a physical site to another but in each step is stored 
in the memories. These are the only quantum systems 
required to possess full coherence robustness along the 
computation. Thus, a number of memories equal to the 
logical inputs suffices for universal computation. 

The present approach is inspired and in turn sheds 
light on other recent proposals for hybrid quantum com- 
puting. For the case of two-dimensional systems, our 
formalism recovers the ancilla-driven model [12 . This 
model, originally derived within a different framework, 
also satisfies properties PI — P3. The main advantage 
here is that, while the relation between the ancilla-driven 
approach and 1WQC is not straightforward [T^EI], our 
derivation immediately enlightens this equivalence. This 
enables the direct application of topological fault-tolerant 
techniques already developed for standard 1WQC to our 
approach [14] , For the CV case, in turn, our formal- 
ism recovers some aspects of the experimental proposal 
of Ref. [15] . There, an implementation of the CV 1WQC 
model was put forward. Nevertheless, this proposal does 
not satisfy property (PI) and requires pulse synchroniza- 
tion, long fiber loops, and non-linear photon-photon in- 
teractions, which impose serious technical obstacles. Our 
proposal alleviates these experimental requirements and 
can be realized using experimentally accessible quantum 
interfaces (16) . In addition, the formalism presented here 
sets a framework within which the proposals of Refs. 
[T2l IT51 [T5] can be understood, and extended to DV sys- 
tems of arbitrary dimension, in a unified way. 

The paper is organized as follows. In Sec. |ll]we re- 
view some basic properties of graph states and their lo- 
cal complementation, for DV and CV systems. There, 
in particular, we also derive an equivalence between the 
SWAP operation and successive local complementations 
for the case of vertices with a single neighbor. In Sec. 
|III| we describe our unified picture for sequential MBQC. 
Physical setups for the implementations of our scheme 



are listed in Sec. |IV| Finally, we present our conclusions 
in Sec. H 



II. GRAPH STATES 

As said, we focus throughout on 1WQC. Its resource 
states are the celebrated cluster states [H 03 Ell HZ], which 
are graph states [7] corresponding to square graphs. Be- 
fore introducing graph states, let us recall the definition 
of mathematical graphs: Consider the union G(y ^ = 
{V, £} of a set V of | V| — N vertices and a set £ of edges 
{j, k} with j, k G V. This can be univocally represented 
by an N x N adjacency matrix T, of elements Tjk £ M^o if 
{j, k} € £ , and Tjk ~ otherwise. If Tjk = 1 V {j, k] € £ 
we call G(v,£) an unweighted graph -or simply a graph- 
arid denote it by G for short; otherwise we write explicitly 
Grp\ and say that the graph is weighted, understanding 
that each edge {j, k} has an associated weight Tjk- 

In both the DV and CV regimes, a graph state \G) 
can be defined as follows: (i) To each vertex j £ V as- 
sociate a quantum system initially in a state given 
by the uniform superposition of all computational states. 
For example, in the case of qubits \+}j = Hj\0}j = 
^(|0)j + \T)j) is the eigenstate of the Xj Pauli operator 

with positive eigenvalue, H is the Hadamard transform, 
and the computational states |0)j and are respec- 
tively the eigenstates with positive and negative eigen- 
values of the Zj Pauli operator. For CV systems in 
turn, \+)j = \0) P] = Fj\0) qj = ^= J R du\u) qj repre- 
sents the zero-momentum eigenstate of the momentum 
quadrature operator pj. In the last, the computational 
states \u) qj are the eigenstates of the position quadrature 
operator <jj, with [qj,pj] — il (we take h = 1 through- 
out), and Fj is the Fourier transform on the j'-th quan- 
tum mode (qumode). (ii) For every edge {j, fc} g £ ap- 
ply a (maximally entangling) controlled- Z gate CZjk to 
registers j and k. For the case of qubits, it is CZjk = 
e tf (i-z ? )®(i-Zfc) < ^ n j s (j e g n j^j on a i so extends directly to 
graph states on d-dimensional systems (qudits) [5]. For 
qumodes in turn, it is CZjk = e %qs ® qk . In all cases it is 

legitimate to write \G) = ri{j,fe}e£ CZ jk <8>jev Fi- 
nally, weighted graph states are obtained by introducing 
multiplicative factors Tij in the exponents of the entan- 
gling operations. 

A. Local complementations of graphs and the 
SWAP operation 

A key family of graph transformations is that of local 
complementations [5J. The definition that we use follows 
closely (and is a simplified version of) the one in Ref. 
[S], more naturally-motivated for weighted graphs. We 
define the local complementation LCn g-\ of Grp), at node 
I G V, and of real weight 5, as the map Grr) — > G(r'n s )■ 
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FIG. 2: Two sequential local complementations: first at ver- 
tex m with weight 1 (left) and then (on the resulting graph) 
at vertex r with weight -1 (right). Black and gray lines corre- 
spond to graph edges with weight 1 and 2, respectively. When 
the initial graph is unweighted and m is the only neighbor 
of r, such composition of local complementations is always 
equivalent to a SWAP operation between vertices m and r. 



Here > T \i,S) jk = r jfe + 5 > if 0> fc } 6 M, and r' m , k = 
Tjk otherwise. The set Ni is the neighborhood of I (in 
G(r)), composed of all vertices n G V connected to I by 
some edge {l,n} € E. In Fig. [2] we show an example 
of two sequential local complementations. An immediate 
but very powerful observation follows for all unweighted 
graphs where one of their vertices -say r— has only one 
neighbor -say to-: The operation LC( m ^ followed by 
LCt r _x) is always equivalent to the exchange (SWAP) 
to -H- r. 

Now, every mathematical transformation LCngs has 
a unitary representation Ulc\i s) (mapping |G(r)) into 
|G(r' (i s ))) which turns out to be local for the aforemen- 
tioned SWAP case. Therefore, if we consider the local 
complementations in Fig. [2] the total transformation, 
for the qubit case, can be expressed in terms of two sim- 
ple local operations: 

U T = Ur®Um= H r ®H m . (1) 

For qumodes, the above transformation is implemented 
similarly as 

U T = U r ® U m = F$ <» F m . (2) 

All these considerations hold true not only for qubits 
and qumodes but also extend straightforwardly to qu- 
dit graph states jS], where the role of the Hadamard 
and Fourier transform gates is played by the discrete 
Fourier transform (recall that the Hadamard is the dis- 
crete Fourier transform for dimension 2). In what follows 
we give explicit proofs just for the qubit and qumode 
cases, i.e. identities ([I]) and 

As we said before, every mathematical transforma- 
tion LC(i t s) has a unitary representation Ulc {is) map- 
ping |G (r) ) into \G [T ' {15) )) = U LC{ls) \G {r) ). Remark- 
ably, this unitary operation turns out to be local. In 
the case of graph-states of qubits for instance, it is 
Ulc u , ±1) = e=Pf^ i£jv -. e^f ^ 7J, which belongs to 
the group of local Clifford operations. Thus, the consid- 
ered succession of local complementations is represented 



by 

= (e**-*- (g) e-*f ) ( e-***- (g) e^ z \ 

= (H r ® H m } (eT^ « e***") . (3) 

Notice that the primed neighborhood N' r refers to the 
neighborhood of r in the graph G(pj ^ ) resulting from 
the first complementation. In addition, the operator 
e - l i x r e ^iz m s i a biii zes a n unweighted graphs having 
a vertex r with a single neighbor to. This can be 
easily checked by noticing that e~ 4 * Xr e l * Zm CZ rm — 
CZ rm e ll i Zm e~' l "^ XrZm , and that the register always 
starts in the state |+). Therefore, restricted to the states 
under consideration, the total transformation ^ is iden- 
tical to 0. 

For qumodes in turn, the corresponding unitary trans- 
formation is Ulc u±1) = e ± 5Pj 1 ® ;eA f. e T 5^ jS], which 
belongs to the group of local Gaussian operations. The 
total composition is now 

U T = e -$# e -39r (g) e i9m e fPm 

= ® F m ) (e*# e~*&) . (4) 

As in the previous case, the operator at the right of the 
Fourier transforms stabilizes the graph states in ques- 
tion, since e^^e~^^CZ rm = CZ rm e^^e l ^ r ^ m , and the 
register is always in the state \p = 0). Thus the total 
transformation Q is identical to 

III. SEQUENTIAL MBQC 

Before introducing our scheme, let us briefly recall the 
conventional approach for the simplest building block of 
the 1WQC, the ID cluster state, also called quantum wire 
[see Fig. |3ja)]: Once the cluster is available the compu- 
tation proceeds by adaptive local measurements along it, 
from left to right in the figure. The measurement-basis 
choices and the length of the wire depend on the par- 
ticular computation to perform. Information flows from 
left to right until all but the particle at the rightmost 
end is measured. The latter encodes the output, and any 
single-particle unitary operation can be implemented in 
such a way. 

As said, the scenario we are interested in involves in- 
stead a static long-lived memory m that sequentially in- 
teracts with flying registers fj. The latter are not re- 
quired to possess robust coherence properties, as they are 
measured immediately after the interaction. All systems 
arc initially in state |+) and each interaction drives & CZ 
gate, so that a graph state among the memory and the 
passing registers is created after each interaction. Our 
aim is to reproduce the conventional 1WQC approach 
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FIG. 3: 1WQC along a quantum wire. The black circles rep- 
resent systems prepared in |+) and the black segments entan- 
gling gates. The node to be measured is depicted in white, and 
the gray circles represent measured nodes, (a) In the original 
scheme measurements are performed from left to right along 
a long linear cluster-state prepared in advance. (6) In the 
considered scenario a static quantum memory m interacts se- 
quentially with the flying register n, then with r2, etc., and a 
star-graph state with m in the center is created (left). Local 
complementations on m and ri swap the positions of m and 
T2 in the associated graph, thus transforming it into a lin- 
ear graph with m at the rightmost end (right), (c) Since all 
the operations (measurements and local complementations) 
involving different flying registers commute, the same com- 
putation can be implemented sequentially by iteration of the 
cycle shown. Notice further that there is no need to actively 
apply U r , since it can always be absorbed in the measurement 
by means of a passive basis-redefinition. 



but with the following variations: only the smallest nec- 
essary pieces of the resource state are created at each step 
of the sequence; and information is stored and processed 
always in the same system, the memory. This would be 
impossible if the flying registers were simply measured as 
they come out from the interaction region. The reason is 
that the considered sequence creates a graph state which 
is not linear but star-like (with m at the center), thus not 
useful for arbitrary single-particle unitary operations. 

A natural way to circumvent this issue is to swap the 
positions (in the underlying graph) of m and before 
the i-th interaction. As said, this can be achieved by 
local complementations as illustrated in Fig. [3^b). For 
the sake of simplicity, let us momentarily assume that 
the measurements are performed only at the end of the 
process. Then, after the interaction between m and the 
first two registers, ri and ra, the graph at the left of 
the figure is created. Before the interaction with the 
next register r^, local complementations at m and at T2 
take the system to the graph state at the right of the 
figure. The process is then repeated with r% and iterated 



throughout, such that a linear cluster with m always at 
the rightmost end is produced. In this way, the ordering 
of measurements in the entire sequence will coincide with 
the flow of information, as in Fig. [3]ja). 

The last ingredient to add for our aim is the obser- 
vation that the operations involving any flying register 
commutes with the ones involving the other flying regis- 
ters. This implies that the reasoning above works equally 
well when measurements are not performed only at the 
end. In particular, it leads us to the sequential imple- 
mentation through the following cycle, schematized in 
Fig. [3^c) : (i) CZ mr . is applied to m and r^. (ii) By lo- 
cal complementations at r; and m their positions in the 
graph are swapped. Finally (Hi) ri is measured, leaving 
m with the quantum information processed exactly as in 
the conventional 1WQC. The process then starts again 
with the next flying register to pass by. 

Universality in 1WQC is achieved with 2D cluster 
states as a resource, which are usually represented by 
square-lattice cluster states. In such a case, horizontal 
connections corresponds to quantum wires, while vertical 
ones to two-body operations [2 . In actual applications, 
however, a fully connected square lattice is not required 
and unnecessary nodes are removed by proper measure- 
ments |2j. Therefore, one can see that the basic building 
block required for sequential MBQC in 2D is represented 
by the graph state at the right of Fig. [4] As shown in 
the figure, each of the graph states from the left hand 
side is related by a local operation (as we have shown 
one can choose the Fourier transform) to the graph state 
at the right. Hence, the sequence is similar to the one 
outlined for the quantum wire: (i) generate one of the 
graph states at the left hand side of Fig. |4j (ii) apply 
the local Fourier transformations, (Hi) implement the lo- 
cal measurements to the registers. Combining these steps 
with the sequential implementation of quantum wires de- 
scribed above, universal quantum computation with n 
logical quantum inputs can be achieved sequentially us- 
ing n quantum memories and n quantum registers that 
are locally measured in each instance of the computation. 
As already mentioned, the registers do not interact with 
each other, but they can act as quantum buses mediating 
the interaction between the memories (as it is the case in 
the bottom- left hand side graph of Fig. [4]). 

To end up with, the method extends to a three- 
dimensional-cluster computation implemented with a 
bidimcnsional cluster of memories. There, since the mea- 
surements are done sequentially layer by layer, the topo- 
logical fault-tolerant techniques [J3] for 1WQC can be 
directly applied in the qubit case. 



IV. IMPLEMENTATIONS 

Several physical setups lie naturally within the range 
of applicability of the ideas presented. For example, 
cavity quantum-electrodynamics (QED) atom-photon in- 
terfaces [TBI US] are perfectly-suited for our scheme in 
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FIG. 4: Sequential MBQC in 2D can be carried out assuming 
that it is possible to build one of the graph states from the 
left side. Applying a Fourier-like transform over each of the 
systems (top) and over rri2 and 62 (bottom) , both graph states 
transform to the one at the right. Then, the computation 
proceeds by measuring n and ri. In the first case there is 
a direct interaction between the quantum memories, while in 
the second one the register, ri, acts as a quantum bus. 



its qubit version: Zero- and one-photon states of fields 
inside microwave resonators can provide the memories, 
and atoms crossing through can play the role of the fly- 
ing registers [HI [20]. In a dual way, ground states of 
static atoms inside optical cavities can be used to en- 
code the memories, and single photons emitted from 
them the flying registers [19H21] , Also artificial-atom- 
photon interfaces can be used analogously, as for in- 
stance quantum-dot [55] and circuit-QED [53] architec- 
tures. Collective atom-light interfaces [TB] in turn give 
yet a further prominent example: Atomic ensembles in 
squeezed spin-coherent states of long-lived ground states 



are sequentially shined through by squeezed-light pulses 
|16) , making thus an excellent candidate for implementa- 
tions in the CV regime. 



V. CONCLUSIONS 

We have introduced a systematic procedure that im- 
plements universal MBQC sequentially. The resulting 
model is less demanding from an experimental point of 
view and maintains the main advantages of the conven- 
tional approach: on-line adaptiveness of entangling oper- 
ations is unnecessary and only local measurements must 
be adapted along the computation. The scheme relies 
mainly on the sequential preparation and manipulation 
of a cluster state on an n-dimensional lattice that, by it- 
erative interactions with auxiliary quantum systems via 
a fixed entangling operation, can perform MBQC cor- 
responding to (n + l)-dimensional lattices. We expect 
that the present proposal contributes to the experimen- 
tal quest for quantum computing and establishes new 
useful grounds for the application of quantum interfaces 
in general. 
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